Breakdown of Scale Invariance in the Phase Ordering of Fractal Clusters 
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Numerical simulations with the Cahn-Hilliard equation show that coarsening of fractal clusters 
(FCs) is not a scale-invariant process. On the other hand, a typical coarsening length scale and 
interfacial area of the FC exhibit power laws in time, while the mass fractal dimension remains 
invariant. The initial value of the lower cutoff is a relevant length scale. A sharp-interface model 
is formulated that can follow the whole dynamics of a diffusion controlled growth, coarsening, frag- 
mentation and approach to equilibrium in a system with conserved order parameter. 
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Non-equilibrium driven dissipative systems relax to 
equilibrium after the driving agent is "switched off' or 
depleted. In complex systems the relaxation dynamics 
can be quite complicated, and it is natural to seek for 
dynamic scaling and universality. An instructive, exactly 
solvable nonlinear example of dynamic scaling in relax- 
ation (coarsening) of rough (self-affine fractal) surfaces 
with nonconservative dynamics is given by the determin- 
istic (undriven) KPZ-equation A much older exam- 
ple is decay of homogeneous and isotropic hydrodynamic 
turbulence [^jQ. Finally, there is an important class of 
relaxation problems related to phase ordering dynamics, 
non-conserved and conserved, in the bulk and on the sur- 
face §§. 

If the system obeys a conservation law, "switching off' 
of the driving agent occurs naturally. There are many 
important non-equilibrium systems that exhibit morpho- 
logical instabilities and ramified growth at an early stage 
of the dynamics, show phase ordering at an intermedi- 
ate stage, and finally approach a simple equilibrium. A 
canonical example is provided by diffusion controlled sys- 
tems, such as deposition of solute from a supersaturated 
solution and solidification from an overcooled liquid. The 
stage of morphological instability and its implications 
have been under extensive investigation [p[-^2[. If some 
noise is present, a fractal cluster (FC) can develop at this 
stage |l2| . The subsequent surface tension driven coars- 
ening of this FC is unavoidable in a closed geometry with 
a finite amount of mass or heat. This stage has not re- 
ceived much attention, with the exception of the paper 
by Irisawa et al jll| where two-dimensional Monte-Carlo 
simulations were performed, and a power law found for 
the perimeter of a DLA cluster versus time. 

We are aware of two additional physical systems with a 
conservation law, for which numerical simulations showed 
nontrivial fractal coarsening dynamics and power laws 



for the cluster perimeter: interface controlled JL3 14 
and surface diffusion controlled Jl3|,[l5| systems. Besides, 
Stokes flow controlled coarsening has been discussed in 
the context of sintering of fractal matter [jl5| . 

No theory is available for any of these fractal coarsen- 
ing systems, except for the very late post-fragmentation 
stage Jig] . On the other hand, a FC is a particular 
case of disordered media with long-range (power-law) 
spatial correlations ]l7|] . The scaling hypothesis (SH) 
(the cornerstone of modern theory of phase ordering [[l8|) 
does not exclude FCs when dealing with long-range cor- 
relations in the initial condition Therefore, one is 
tempted to employ the SH and calculate the growth ex- 
ponents for the coarsening of FCs. We start with these 
simple calculations. Then we report our simulations of 
the diffusion controlled coarsening of a DLA aggregate, 
as described by the Cahn-Hilliard (CH) equation. Hav- 
ing measured, for the first time, the dynamics of the pair 
correlation function (which is very close to the average 
mass density) of the FC, we show that the SH is invalid. 
On the other hand, we find that a characteristic coars- 
ening length scale and interfacial area of the FC exhibit 
power law dynamics (with a new growth exponent), while 
the fractal dimension remains invariant (on an interval 
of scales shrinking with time). The initial value of the 
lower cutoff of the FC is shown to be an additional rele- 
vant length scale. Finally, a minimalistic sharp-interface 
model is presented that can follow the whole dynamics 
of the diffusion controlled system: an unstable growth, 
coarsening, fragmentation and approach to equilibrium. 

Let the initial state of a conserved system represent 
a single-connected, statistically homogeneous self-similar 
mass fractal of the minority phase, characterized by the 
fractal dimension D on an interval of scales between the 
lower cutoff Iq and upper cutoff L. We start with a simple 
coarsening scenario [pl|pO| that is required by the SH. It 
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assumes that the fractal dimension of the coarsening clus- 
ter remains constant on a (shrinking) interval of scales 
between the time-dependent lower and upper cutoffs, l(t) 
and L{t). The interfacial area A and total mass M of the 
FC are estimated as [fL7f 



A~I d - 1 (£/0 D and M~l d (L/lY 



(1) 



respectively, where d is the embedding Euclidean dimen- 
sion. Mass conservation yields L ~ l( D ~ d )/ D [^FJ. Now 
assume that l(t) -t 1 / 2 . Then we find the following scal- 
ing laws: A(t) ~ f- 1 /* and L(t) ~ t ( -D~ d )/ Dz . The scal- 
ing of Lit) describes shrinking of the FC in the process 
of coarsening |Q . 

Already at this stage a discrepancy appears: no shrink- 
ing has been observed in any direct numerical simula- 
tions of coarsening of FCs |l3|-[l5|]. This gives a strong 
evidence for breakdown of scale invariance pi]] . On the 
other hand, power laws for A(t) reported in Rcfs. 13 lj| 
indicate that the problem might possess scaling behavior 
of a more complicated nature. 

To clarify the matter, we performed more detailed nu- 
merical simulations of a diffusion-controlled system. In 
addition to L(t) and Ait), we followed the evolution of the 
equal-time pair correlation function (which is very close 
to the average mass density of the FC, so we will not dis- 
tinguish between them). Having measured it, one could 
find the mass fractal dimension and coarsening length 
scale of the FC for every moment of time. 

If one remains, for one more moment, within the frame- 
work of the SH, one can easily predict the dynamics of 
the mass density p(r,t). At distances r <C l(t) from a 
(typical) reference point inside the cluster, the cluster 
is non-fractal: p(r, t) ~ const. At distances interme- 
diate between I and L p(r, t) = a(t) r D ~ d , where a(i) 
is a function of time. Matching these two asymptotics, 
we have ait)l D ~ d = const and hence a(t) ~ t( d ~ D ^ z . 
Therefore, for l(t) < r < L(t) the SH predicts p(r,t) ~ 
[r /t 1 / z ) D ~ d , a simple self-similar expression. It is the ab- 
sence of this self-similarity that will enable us to utterly 
disprove the SH. 

We concentrated on the diffusion controlled coarsen- 
ing and employed the CH equation, a standard model of 
phase ordering with a conserved order parameter (COP) 
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V 2 (V 2 m + u-u 3 ) 



= 0. 



(2) 



Eq. (g) was discretized and solved on the domain < 
x < 512 , < y < 512 with periodic boundary condi- 
tions. We used an explicit Euler integration scheme to 
advance the solution in time, and second order central 
differences to discretize the Laplace operator. With a 
mesh size Ax = Ay = 1 no preferred directions emerged 
in the computational grid, due to the truncation errors; a 
time step At = 0.05 was required for numerical stability. 



The accuracy was monitored by checking the mass con- 
servation that was verified in all the simulations within 
0.01 %. 

We chose a DLA cluster |22j as the initial condition. 
The fractal properties of DLA clusters are somewhat 
more complex that those of a simple self-similar frac- 
tal pa]. However, it is a DLA- like FC that can develop 



during the diffusion controlled growth |12|, so this choice 
is physically motivated. The initial clusters (like the one 
shown in Fig. la), with radius of order 250, were pre- 
pared by a standard random-walk algorithm on a two- 
dimensional square grid. To prevent fragmentation at 
an early stage of the coarsening process, we followed the 
technique of Irisawa et al. Jl^,|lJ]: the aggregates were 
thickened by an addition of peripheral sites. The mass 
fractal dimension was determined from the mass-radius 
relation E^] and ranged from 1.67 to 1.72. 

We identified the cluster as the locus where u(r, t) > 
0. The coarsening process was followed up to a time 
t = 5, 000. Typical snapshots of the coarsening pro- 
cess are shown in Fig. 1. One can see that smaller fea- 
tures of the FC are "consumed" by larger features, while 
the global structure of the cluster is not affected. To 
characterize the coarsening process, the following quan- 
tities were sampled and averaged over 10 initial con- 
figurations: (1) the gyration radius of the cluster, (2) 
the circularly averaged pair correlation function <?(r, t) — 
([u(r' , t) + 1] [u(r' + r, t) + 1]), (3) the cluster perimeter 
Ai(t), defined as the sum of | Vti(r, t) | 2 over the whole 
domain, and (4) the cluster perimeter ^(t), defined as 
the number of broken bonds between the aggregate sites. 

The gyration radius of the FC has been found to re- 
main constant within possible logarithmic corrections. 
Evolution of g(r, t) is shown in Fig. 2. One can see that 
coarsening affects only the smallest lengths, while the 
intermediate-distance power-law part remains "frozen" . 
It is evident that g(r, t) does not acquire a self-similar 
form, so there is no scale invariance. On the other hand, 
the mass fractal dimension remains invariant on an in- 
terval l c {t) <C r <C L. The dynamics of the coarsening 
length scale l c (t), extracted for each moment of time from 
the slope of the linear part of g(r, t) (the Porod law M), 
are shown in Fig. 3. The late-time behavior of the slope 
versus time shows a power law: t~ ai with a\ = 0.19. 
Therefore, l c {t) ~ t" 1 , and the corresponding growth ex- 
ponent z\ = 1/ai is close to 5 (and not to 3 as could be 
expected for a diffusion controlled system with a COP 
S). Fig. 4 shows the dynamics of the cluster perimeter 
estimates A\it) and ^(t). The long-time dynamics of 
each of them is describable by a power law t~ a2 , with 
a-2 = 0.19 for Ai, and 0.20 for A^. The corresponding 
result of Monte-Carlo simulations [^3) was slightly differ- 
ent: 0.22 — 0.24. The close proximity of the exponents 
o.\ and ai gives evidence that it is a single exponent. 
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FIG. 1. Evolution of a DLA cluster undergoing coarsening 
in a conserved, diffusion controlled system. The upper row 
corresponds to t — (left) and 34.7 (right), the lower row to 
t = 329.3 (left) and 4,900 (right). 
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FIG. 2. Dynamics of the pair correlation function g(r,t) 
for time moments t = (solid line), 516.5 (dotted line) and 
4,900 (dashed line). 

Absence of scale invariance means presence of an addi- 
tional length scale. Function g(r, t) gives evidence for the 
nature of this length scale. For our mass fractal at t = 
we have g(r, t = 0) <~ (r/lo) D ~ d in the fractal region 
la <C r « i. Preservation of the power-law part of g with 
time (Fig. 2) implies that the same asymptotics holds, on 
a shrinking interval of radii, for t > (until fragmenta- 
tion). That is, the small intrinsic length scale la remains 
relevant. How does it show up in the phenomenology of 



coarsening? Fig. 1 gives evidence that (i) the FC can be 
regarded as a set of "bars" , and (ii) the characteristic bar 
length lb grows in time faster than the bar width (iden- 
tified with l c ). The area of a single bar should scale like 
lb l c , hence the total area of the FC is lb l c (L/lb) D . This 
quantity must be equal to the initial value of the FC area, 
ll {L/h) D . This yields l b ~ l (l c /h) 1,{ - D - 1) ~ t^/^K 
We will finish this Letter with formulating a sharp- 
interface model that can describe the whole diffusion- 
controlled dynamics, from the stage of growth through 
coarsening and fragmentation to the final equilibrium. 
Consider a number of (possibly multiple-connected) mass 
clusters characterized by a set of their (moving) interfaces 
7i. Let now u(r,t) be the mass concentration of the so- 
lution normalized to the (constant) density of solute in 
the compact solid phase. The field u in the liquid phase 
is governed by the diffusion equation 
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FIG. 3. The slope of the linear part of g(r,t) versus time, 
and its power-law regression. 
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FIG. 4. Two estimates for the FC perimeter versus time: 
Ai (triangles) and A2 (circles). 



du ~ 



(3) 



in a finite c?-dimensional domain. We specify a no-flux 
boundary condition, V„it |r= on the external bound- 
ary r, where index n stands for the normal component 
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of a vector. Assuming that each of the interfaces 7, is in 
local thermodynamic equilibrium, we employ the Gibbs- 
Thomson relation u | 7( = uo(l + Ao«i), where uq is the 
(normalized) equilibrium concentration of the solution in 
the bulk, Ao is the capillary length and k.; is the local 
curvature for d = 2, or the mean curvature for d > 2. 
(We limit ourselves to an isotropic surface tension.) Fi- 
nally, mass conservation at each of the moving interfaces 
yields the well-known relation for the normal speed: 



1-U 



(4) 



It is easy to check that this model preserves the total 
mass of the solute. In the normalized form 



u dr = const , 



(5) 



where ft c is the total volume (area) of the solid phase, 
while Q denotes the region unoccupied by the solid phase. 
This important conservation law does not appear in the 
more traditional theoretical formulations of the diffusion- 
controlled growth problem ]8|-[l2]|, where an "infinite" 
system is studied, and the boundary condition corre- 
sponding to a constant (positive) flux or constant su- 
persaturation at r — > 00 is used. Notice that, even in 
the limit of strong diffusion, it is the full diffusion equa- 
tion (rather than its Laplace's equation limit) and no-flux 
condition on T that provide the conservation law. Also, 
the usually small term u in the denominator of Eq. (^) 
should be kept to get Eq. (||) right. 

In summary, we have demonstrated that diffusion con- 
trolled phase ordering of FCs is not a scale-invariant pro- 
cess. In spite of this, the problem possesses non-trivial 
scaling properties: the coarsening length scale and inter- 
facial area of the FC exhibit power laws in time (with a 
new growth exponent), and the mass fractal dimension 
remains invariant. An additional small intrinsic length 
scale (the initial value of the lower cutoff) remains rele- 
vant until the fragmentation stage. We believe that these 
findings apply to other coarsening mechanisms as well. 
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